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1.  INTRODUCTION 


Traditionally,  quality  control  charts  have  been  designed  with  re¬ 
spect  to  statistical  criteria  only,  and  the  control  methodology  is  based 
on  the  independence  and  normality  of  serial  samples.  At  first  the  pro¬ 
duction  process  is  assumed  to  be  characterized  by  a  single  in-control 
state.  For  example,  if  the  process  has  one  measurable  quality  character¬ 
istic,  then  the  in-control  state  will  correspond  to  the  mean  of  this 
quality  characteristic  when  no  assignable  cause  is  present.  ,  ^  -  1 


Now  we  consider  the  model: 


xt  ■ 11  +  5t 


(l.D 


where  p  is  a  constant,  is  an  error  and  is  t-th  observation.  It  is 

of  interest  to  select  the  sample  size,  statistical  characteristic  of 

and  control  limits  so  that  the  power  of  the  test  to  detect  a  particular 

shift  in  the  quality  characteristic  and  the  type  I  error  probability  are 

equal  to  specified  values.  Usually,  t  =  1,2,  ...  is  considered  to 

be  independently  normally  distributed  with  zero  mean  and  common  variance 
2  2 

a  ,  where  a  is  known  or  unknown.  In  this  case,  consideration  of  statis¬ 
tical  criteria  and  practical  experience  have  led  to  general  guidelines 
for  the  design  of  control  charts  resulting  in  widespread  use  of  samples 
of  size  5,  three-sigma  control  limits,  and  a  sampling  frequency  of  one 
hour  for  the  X-chart  (see  Duncan  [3]).  In  the  sequel,  we  set  j  ■" 


*  - 1  l  V 

i  =  l 
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,  -urt er  of  samples.  Hereafter  we  use  X 
;^rerJi  3.  quantity  control  limits  on 
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(1.3) 


“«*  ■  *  1  -  ; j -*  1  >«iation  is  known  or  estimated. 

limits  on  the  standard  deviation  are 
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(1.4) 


(1.5) 


'"on  or  estimated,  where 


r  , 


(1.6) 


literature  {see  Grant  [4]).  The  process  will  be  considered  under  con¬ 
trol  if  the  estimate  of  the  mean  and  the  estimate  of  the  standard  de¬ 
viation  of  the  process  remain  within  prescribed  control  limits  above. 

In  practice,  a  number  of  data  sets  in  economics,  business,  enqineerim 
and  the  natural  science  often  are  present  in  the  form  of  time  series. 

In  other  words,  the  observations  are  dependent,  i.e.,  ?t's  of  model  (1.1) 
are  not  white  noise;  for  example,  U^,  t  =  0,±l,  ...}  is  an  autoregressive 
moving  average  (ARMA)  with  order  (p,q).  So  the  problem  is  how  to  deter¬ 
mine  3a  control  limits.  Stamboulis  [7]  studied  AR(1)  with  parameter  a. 
Vasilopoulos  [8]  extended  Stambouli's  results  to  ARMA(p,q)  model.  Vas¬ 
il  opoul  os  and  Stamboulis  [9]  together  investigated  the  case  AR(2)  = 
ARMA(2,0) .  It  is  different  from  classical  control  factors.  How  diff¬ 
erent  it  is  depends  on  the  stochastic  properties  of  the  process.  Since 
the  method  is  similar,  we  only  discuss  AR(2). 


2.  CONTROL  CHART  ON  THE  SECOND  ORDER  AUTOREGRESSIVE  MODEL 
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In  model  (1.1),  assume  Is  an  AR(2)  model,  that  is 

h  =  al5t-l  +  a2ct-2  +  et  (2J) 

p 

where  {e+}  is  a  white  noise  series  with  Ee.  =  0  and  V(e+)  =  a  ,  a,  and 

l  t  L  C  I 

a2  are  constants.  For  stationarity  of  AR(2),  it  is  necessary  that  the 
roots  of  the  characteristic  equation  of  the  AR(2)  process 

♦  (B)  =  1  -  c^B  -  a2B2  =  0  (2.2) 


must  lie  outside  the  unit  circle,  which  is  equivalent  to  the  alpha  co¬ 
efficients  being  in  the  triangular  region: 

a2  +  °1  <  1  *  a2  ”  al  <  1  *  -1  <  a2  <  1 


(see  Box  and  Jenkins  [2]).  The  variance  of  the  AR(2)  process  is  given 
by 


2 
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I(l-a2)2  -  a2] 


(2.3) 


Suppose  ,  k  =  0,±l,  ...,  are  the  autocovariance  functions  of  the 

2  2  — 

AR(2)  process,  then  a  and  the  variance  a—  of  the  sample  mean  x  are 

given  in  terms  of  y ^  by 

O2  =  Ynl  (2.4) 


4*1Cy0*2J,  O'nK1-  (2-5) 

In  order  to  evaluate  the  control  limits  for  Y,  we  need  to 

O 

evaluate  at.  .  This  can  be  accomplished  by  (2. 3)- (2. 5)  if  ,  a2  and 

2  2 
a  are  known.  Therefore, if  the  process  variance  a  is  known,  the  control 

e 


limits  of  the  model  described  by  (1.1)  and  (2.1)  is  modified  to 


X  ±  A(a-|  ,a2 ,n)a 


(2.6) 


where 

A(a-j  ,<*2 >n)  =  X1>/2(a1  ,a2»n)  ~  ,  (2.7) 

n-1  . 

X(a1,a2,n)  =  1+2  l  (l-£)bt,  (2.8) 

t "  1 

bt  =  VY0*  (2o9) 

o 

The  expression  of  x{a.|  ,a2 ,n)  is  based  on  the  expression  of  a-  . 

In  order  to  construct  the  7-chart  when  the  standard  deviation  is 
unknown,  we  must  evaluate  the  auxiliary  parameter  C2(ct.j  ,a2,n)  first, 

which  is  also  needed  to  construct  control  limits  of  the  standard  devia- 

2  2  2 
tion.  Since  S  /ES  is  distributed  as  xn»  we  get 


E(S) 

Hence, 


fj  r(y)  (  ?  n-1  .  \l/2 

C2 (otl  ,a2,n)q  =  /n  ^Twiyl1  "  ~  n^t)  a>  (2J0) 


/  o  n-l  t  \1/2 
C2(0l,a2,n)  =  C2-^l -CT^(1 --)btJ  . 


(2.11) 


To  obtain  an  approximate  expression  for  E(S),  we  use  this  ex¬ 
pression  of  E(S): 


E(S)  =  zJ?  =  -  -^4)  - 

V  8{ES2)2/ 


4? 

sT^TT 


(2.12) 


In  general,  the  last  term  of  (2.12)  is  smaller  than  the  first  term. 
For  example,  let  n  =  5.  Then  the  last  term  of  (2.11)  is  0.002144* 


so  we  can  omit  this  term  and  get  the  approximate  formula: 


E(S)  *  ^  . 

\  8(ES2)2 


The  expected  value  and  variance  of  S  can  be  obtained  by 


ES2  = 


/  X(apa2*n)\ 

v  "  / 


and 


O  1  ^  n  n 

Var(S2)  =  X  I  v.t+2 
n  t=l  x= 1 


,  n  n 

T  I  I  Yt_ 

n  t=l  t=1  r  T 


-i  2 


(2.13) 


(2.14) 


‘  ^  J,  L  Yt-xYt-Y‘ 

n  t=  I  t- I  y  l 


(2.15) 


2 

But,  from  the  expression  of  Var($  ),  the  complexity  of  (2.13)  is  not 
better  than  (2.10). 

By  replacing  A-j  (a.j  ,a2,n) ,  C^a-i  ,ct2>n) ,  B^  (a^,a2,n)  for  A-j,  and 
B^ ,  i  =  l,2,3,4  of  (1.3)  and  (1.6)-(1.8),  respectively,  the  modifica¬ 
tion  of  control  chart  limits  in  an  AR(2)  model  is  obtained. 

The  substantial  ranges  in  the  values  of  x(a-|,a2,n)  and  C2(a^,a2,n) 
greatly  affect  the  control  factors.  Vasilopoulos  and  Stamboulis  [9] 
gave  an  example  to  illustrate  this  result. 


3.  CONTROL  CHARTS  IN  MULTIVARIATE  CASE 


Now  we  consider  multivariate  case.  The  model  in  this  case  is 


*t  = 


+  St 


(3.1) 


where  y  is  a  mxl  constant  vector,  {5  ,  t  =1,2,  . . .  >  is  a  m-dimensional 
stationary  process  with  zero  mean  vector.  Set 


1  - 1 


nt£,Sl*  s  =  srr  .Hii-iKSi-x)', 


where  the  prime  means  transpose  of  a  matrix  (or  vector).  Let  x  be  the 
global  mean  over  several  subgroups  of  size  n,  and  ST  be  the  pooled  sample 
covariance.  It  is  well-known  that  if  Ut)»  t  =  1,2,  ...  is  a  series  con¬ 
sisting  of  white  noise  with  distribution  N  (o,A),  the  X-control  chart 

m  ~  - 

has  been  studied  by  many  authors  (for  example,  Ghare  and  Torgrersen, 
Jackson  and  et  al).  From  the  facts  that  (x-  y) 1  (^-a)-1  (x-  y)  and 
(x-y)'(-S)  (x-y)  are  distributed  as  xm  and  Hotelling  T  -statistics 

respectively,  we  can  construct  the  quality  control  region  on  x  based  on 
A  is  known  or  unknown.  When  a  is  known,  the  control  region  is 


D  =  f x :  n(x  -  x) 'A-1  (x  -  x)  <  x^(o)}. 


(3.2) 


This  is  an  elliptical  region.  When  A  is  unknown,  the  control 
region  is 


D  = 


l*:  ^(x  -  *)  'S'1  (X  -  X)  lFm>n.m+1(a)}. 


(3.3) 


But  in  practice,  { }  is  not  generally  a  white  noise  series.  When 
{ )  is  serially  dependent  and  described  by  p-dimensional  ARMA(p,q) 
model,  the  control  region  of  mean  vector  will  be  modified.  For  simplicity. 


8 

we  will  discuss  p-dimensional  AR(P)  model.  Let 

Xt  =  y  +  Ct,  (3.4) 

St  =  BlSt-l  +  +  BpCt-p  +  if  (3’5) 

where  B. ,  i  °1 . p,  are  mxm  matrices  and  a  white  noise  series  with 

distribution  N^CO.E),  £  >  0.  Furthermore , we  can  also  generalize  (3.5)  to 
the  following  form: 

St  ’  Mt-l  +  •••  +  Bp6t-p  *  fl£t.  (3.6) 

where  A:  mxr,  et  -  Nr(o,£)  such  that  AeA'  >  0.  The  model  described  by 

(3.5)  and  (3.6)  is  often  met.  For  example,  in  the  production  of  synthe¬ 
tic  fiber  the  tensile  strength  x-|  and  diameter  x^  may  be  equally  impor¬ 
tant  quality  characteristics.  Their  flucuations  mainly  result  from  mois¬ 


ture,  then,  in  proper  productive  process,  (x^.x^)'  may  be  described  by 

(3.5)  and  (3.6).  Here  we  only  discuss  the  model  (3.5)  because  the  method 
treating  model  (3.6)  is  the  same  as  (3.5). 

It  is  well-known  that  the  necessary  condition  that  the  AR(P)  model 

(3.5)  is  stationary  is  all  the  roots  of  determinant  of  (XPI  -  xp~1B1  - 
...  -  XBp_.j  -  Bp)  lie  within  the  unit  circle.  Set 

Ak  =  ^-t  ~  ^-t+k  ~  ’ 


then,  there  also  exist  "multivariate  Yule  Walker"  equation: 


Hence,  if  B. ,  i  =  1 ,  . .. ,  p,  and  z  are  i denti fied  from  data,  then  al  1 
Ak.  k  =  0,  1 ,  . . . ,  can  be  calculated  from  (3.7)-(3.9).  Furthermore,  the 
covariance  matrix  of  7  can  be  obtained: 


Ax-i(vX(’-X*n>)  • 


(3.10) 


where  Aq  is  the  covariance  matrix  of  x^. 

C  -i  n  ro  {  v  _  n  W  A  ^  ^  v  ,  i  \  ■»  -i  , ,  “  -C  ^  ■ 


Since  (x  -  p) ' A —  (x  -  y)  is  ax  distribution  with  degrees  of  freedom 
we  can  get  the  control  region  of  mean  vector  within  an  elliptical  region: 


D  =  {x:  (x-  x)'Al1(x-  x)  <  x^(a)>. 


(3.11) 


Notice  that  if  §t  =  in  (3.5),  i.e.,  our  process  is  classical, 
then  the  control  region  described  by  (3.11)  is  the  same  as  the  one  de¬ 
scribed  by  (3.2).  When  p  f  0,  these  elliptical  regions  described  by 
(3.2)  and  (3.11)  are  different  from  the  lengths  and  directions  of  their 
major  axes. 
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